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1.  Introduction 


Precision  munitions  research  is  currently  an  area  of  focus  throughout  the  United  States  Armed 
Services.  Within  the  Army,  rolling  munitions  using  movable  lifting  surfaces  for  maneuver 
control  have  been  and  are  continuing  to  be  developed.1,2  However,  the  maneuverability  of  such 
munitions  is  limited  to  the  munition  being  positioned  in  the  proper  orientation  to  move  in  the 
desired  direction.  Research  shows3  that  a  non-rolling  airframe  using  4  independent  canards  in  a 
skid-to-turn  maneuver  is  the  most  effective  method  by  which  to  guide  a  munition  with  movable 
lifting  surfaces.  As  the  armed  services’  desire  to  engage  moving  targets  as  well  as  stationary 
targets  increases,  the  need  for  understanding  the  aerodynamics  associated  with  highly 
maneuverable  airframes  also  increases. 

To  that  end,  the  2013  US  Army  Research  Laboratory  (ARL)  Low  Cost  Hyper  Accurate  Weapon 
(LCHAW)  mission  program  seeks  to  develop  a  highly  maneuverable  projectile  that  is  able  to 
acquire,  track,  and  intercept  a  moving  target  as  the  natural  extension  of  the  Very  Affordable 
Precision  Projectile  (VAPP)  program  that  was  successfully  able  to  guide-to-hit  a  stationary 
target.  This  effort  requires  research  in  a  few  key  areas:  high  maneuverability  airframe,  maneuver 
technology,  guidance  and  flight  control,  and  image-based  navigation,  all  at  low  cost. 

Prior  to  engaging  in  the  research  on  the  high  maneuverability  airframe  and  maneuver  system,  a 
demonstration  platform  needed  to  be  selected.  Choosing  a  demonstration  platform  for  the 
moving  target  mission  program  required  a  projectile  configuration  that  would  meet  stability, 
range,  and  maneuverability  requirements  for  the  demonstration,  while  remaining  simple  enough 
so  that  quality  research  could  be  engaged.  Initially,  the  83-mm  shoulder-launched  munition 
under  development  at  the  Armament  Research,  Development,  and  Engineering  Center 
(ARDEC),  Picatinny  Aresenal,  NJ,  was  chosen.  However,  previous  work4  found  that  that 
munition  would  not  meet  the  requirements  of  the  current  demonstration  platform.  Still,  it  was 
decided  that  this  83-mm  shoulder-launched  munition  would  serve  as  the  starting  point  for  the 
design  because  it  had  many  positive  characteristics. 

The  front  end  of  the  shoulder-launched  munition  was  satisfactory  from  a  flight  perspective  and 
provided  enough  volume  for  guidance  sensors  and  electronics.  However,  the  rear  end  of  the 
munition  did  not  meet  the  aerodynamic  requirements.  Therefore,  it  was  decided  that  the  rear  end 
of  the  shoulder-launched  munition  would  be  modified  to  that  of  a  standard,  artillery-like 
boattail.4  While  the  previous  study4  included  optimization  of  both  the  fins  and  canards,  the 
interaction  effects  between  them  were  not  considered.  Previous  research  has  shown  that  these 
interaction  effects  are  critical  to  the  performance  of  the  projectile  because  the  downwash  of  the 
canard  vortices  can  reduce  the  effectiveness  of  the  fins,  which  could  have  detrimental  effects.2,5  7 


1 


This  report  describes  the  computational  fluid  dynamics  (CFD)  simulations  and  follow-on 
analyses  completed  to  optimize  the  entire  munition’s  configuration,  including  the  maneuver 
system  and  interaction  effects.  The  resulting  configuration  will  be  used  as  the  demonstration 
platform  for  the  remainder  of  the  research  under  the  LCHAW  mission  program. 


2.  Geometry 


Three  body-tail-canard  configurations  were  investigated  during  the  course  of  this  study: 
Configuration  la,  Configuration  lb,  and  Configuration  2.  All  aerodynamic  coefficients  were 
determined  for  simplified  geometries;  all  gaps  and  slots  in  the  body  that  would  be  needed  to 
accommodate  the  fins  and  canards  were  not  included.  Figure  1  shows  a  representative  body-tail- 
canard  configuration.  Each  geometric  model  was  created  in  SolidWorks8  and  has  a  body 
reference  diameter  of  83  mm  (1  caliber).  The  body  geometry,  which  did  not  vary  among 
configurations,  consisted  of  a  hemispherical  nose  cap;  a  nearly  cylindrical  body  320.38  mm  (3.86 
cal.)  long;  and  a  66.4  mm  (0.8  cal)  long,  7°  boattail.  The  overall  length  was  427.45  mm  (5.15 
cal.).  During  the  course  of  the  study,  the  projectile  center-of-gravity  ( Xcg ),  the  number  of  fins, 
the  fin  span,  chord,  and  geometry,  and  the  canard  dither  axis  location  (xc)  and  span  (but  not 
chord)  were  varied  between  the  configurations.  A  summary  of  the  3  configuration  designs  can  be 
found  in  Table  1.  Specifics  of  each  configuration  are  discussed  in  Section  4. 


Fig.  1  Representative  body-fin-canard  geometry  used  during  the  course  of  this  study 
Table  1  Configuration  design  summary 


Configuration 

Xcg  (mm 
From 

Nose) 

No. 

of  Fins 

Fin  Span 
(mm) 

Fin 

Chord 

(mm) 

Fin  Geometry 

xc  (mm 

Forward 

otXc„) 

Canard 

Span 

(mm) 

Canard 

Chord 

(mm) 

la 

200 

4 

188.4 

22.0 

double  bevel 

68.89 

240.48 

18.86 

lb 

200 

4 

188.4 

22.0 

double  bevel 

91.30 

192.72 

18.86 

2 

218 

8 

199.9 

19.05 

bent  fin  tip 

73.79 

192.72 

18.86 
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3.  Solution  Techniques 


3.1  Computational  Aerodynamics 

CFD  simulations  were  completed  on  all  3  geometries  to  characterize  the  aerodynamic  behavior 
of  the  designs  prior  to  wind  tunnel  or  flight  tests.  This  section  briefly  describes  the  methodology 
used  for  the  simulations. 

3.1.1  Solver 

The  double-precision  solver  of  a  commercially  available  code,  CFD++  vl2.1.1,9  is  used  for  the 
CFD  simulations.  The  basic  numerical  framework  in  the  code  contains  unified-grid,  unified- 
physics,  and  unified-computing  features.  A  brief  synopsis  of  this  framework  and  methodology  is 
given  below.  Additional  information  is  available  in  Pulliam  and  Steger10  and  in  Peroomian  et 
al.11 

The  3-dimensional  (3-D)  Reynolds-averaged  Navier-Stokes  (RANS)12  equations  are  solved 
using  the  following  finite  volume  method: 

jJvWd.V  +  $[F-G]  -dA  =  fvHdV,  (1) 

where  W  is  the  vector  of  conservative  variables;  F  and  G  are  the  inviscid  and  viscous  flux 
vectors,  respectively;  H  is  the  vector  of  source  terms;  V  is  the  cell  volume;  and  A  is  the  surface 
area  of  the  cell  face. 

The  numerical  framework  of  CFD++  is  based  on  the  following  general  elements:  1)  unsteady 
compressible  and  incompressible  Navier-Stokes  equations  with  turbulence  modeling  (unified- 
physics);  2)  unification  of  Cartesian,  structured-curvilinear,  and  unstructured  grids,  including 
hybrids  (unified-grid);  3)  unification  of  treatment  of  various  cell  shapes,  including  hexahedral, 
tetrahedral,  and  triangular  prism  (i.e.,  3-D)  cells;  quadrilateral  and  triangular  (i.e.,  2-dimensional 
[2-D])  cells  and  linear  (i.e.,  1-dimensional  [1-D])  elements  (unified-grid);  4)  treatment  of 
multiblock  patched  aligned  (nodally  connected),  patched-nonaligned,  and  overset  grids  (unified- 
grid);  5)  total  variation  diminishing  discretization  based  on  a  new  multidimensional  interpolation 
framework;  6)  Riemann  solvers  to  provide  proper  signal  propagation  physics,  including  versions 
for  preconditioned  forms  of  the  governing  equations;  7)  consistent  and  accurate  discretization  of 
viscous  terms  using  the  same  multidimensional  polynomial  framework;  8)  point-wise  turbulence 
models  that  do  not  require  knowledge  of  distance  to  walls;  9)  versatile  boundary  condition 
implementation  that  includes  a  rich  variety  of  integrated  boundary  condition  types  for  the  various 
sets  of  equations;  and  10)  implementation  on  massively  parallel  computers  based  on  the 
distributed-memory  message-passing  model  using  native  message-passing  libraries  or  Message 
Passing  Interface  (MPI),  Parallel  Virtual  Machine  (PVM),  etc.,  (unified-computing). 
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The  code  has  brought  together  several  ideas  on  convergence  acceleration  to  yield  a  fast  steady- 
state  solution  methodology  for  all  flow  regimes.  The  approach  can  be  labeled  a  “preconditioned- 
implicit-relaxation”  scheme.  It  combines  3  basic  ideas:  implicit  local  time-stepping,  relaxation, 
and  preconditioning.  Preconditioning  the  equations  ideally  equalizes  the  eigenvalues  of  the 
inviscid  flux  Jacobians  and  removes  the  stiffness  arising  from  large  discrepancies  between  the 
flow  and  sound  velocities  at  low  speeds.  The  use  of  an  implicit  scheme  circumvents  the  stringent 
stability  limits  suffered  by  their  explicit  counterparts,  and  successive  relaxation  allows  update  of 
cells  as  information  becomes  available  and,  thus,  aids  convergence.  The  suggested  ramping  of 
the  Courant-Friedrich-Lewy  (CFL)  number  within  CFD++  also  speeds  convergence. 

Only  the  subsonic  regime  (0.3  <  M  <  0.8)  is  considered  in  this  report,  where  the  suggested  CFL 
ramping  is  from  1.0  to  100.0  over  100  iterations.  However,  depending  on  the  angle  of  attack 
being  investigated,  the  ramping  of  the  CFL  number  had  to  be  limited  to  a  value  of  25  or  50. 
Second-order  discretization  was  used  for  the  flow  variables  and  the  turbulent  viscosity  equation. 
Turbulence  closure  is  based  on  topology-parameter-free  formulations.  Only  1  of  the  available 
turbulence  models  within  CFD++  was  used  in  the  current  study,  based  on  the  results  of  the 
turbulence  study  in  Silton  and  Fresconi4:  the  Goldberg  3-equation  k-s-R  turbulence  model.13 
This  3-equation  model  solves  the  transport  equations  for  undamped  eddy  viscosity,  R,  in  addition 
to  turbulent  kinetic  energy  ( k )  and  its  dissipation  rate  (s).  This  solution  methodology  accounts  for 
non-equilibrium  conditions  and  avoids  freestream  turbulence  decay  under  shear-free  flow 
conditions. 

A  tool  exists  within  CFD++  to  recommend  values  for  the  initialization  of  the  turbulence  transport 
based  on  user  specified  free-stream  turbulence  intensity  (or  the  turbulent-to-molecular  viscosity 
ratio)  and  the  turbulence  length  scale,  which  the  model  requires.  For  the  current  problem,  the 
turbulence  intensity  was  set  to  2%  and  the  turbulent-to-molecular  viscosity  ratio  was  set  to  50 
since  the  length  scale  was  not  known. 

3.1.2  Numerical  Grids 

All  grids  used  in  the  computational  aerodynamics  simulations  were  created  using  MIME  v4.1.14 
MIME  is  an  unstructured  mesh  generator  that  allows  triangular  or  quadrilateral  dominant  cells 
for  the  surface  mesh.  Once  an  adequate  surface  mesh  is  generated,  prism  layers  can  be  specified 
and  are  created  when  the  volume  mesh  is  generated.  MIME  saves  a  parameter  file  with  the 
surface  names  and  desired  cell  sizes  that  can  be  utilized  for  creating  new  meshes  with  similar 
geometries  (as  in  the  present  study).  The  user  only  needs  to  import  the  new  geometry,  open  the 
old  parameter  file,  and  assign  the  geometry  as  appropriate.  The  grid  employed  here  is  based  on 
the  parameter  file  used  to  create  the  final  mesh  in  Silton  and  Fresconi.4  As  such,  no  further  mesh 
refinement  study  was  conducted. 

Figure  2  shows  the  limits  of  the  outer  boundary  of  the  mesh.  It  extends  approximately  35  body 
lengths  from  the  projectile  in  all  directions.  Cylindrical  density  boxes  were  placed  in  the  wake  of 
the  projectile  and  the  fins  to  ensure  proper  resolution  in  this  area.  An  additional  density  box  was 
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placed  between  the  canards  and  the  fins  to  ensure  that  the  canard  tip  vortices  would  not  dissipate 
too  much  before  reaching  the  fins.  The  surface  mesh  and  boundary  layer  growth  was  restricted  to 
a  ratio  of  1.2  or  less.  The  growth  ratio  was  relaxed  to  2.0  when  the  remainder  of  the  volume 
mesh  was  created.  The  prism  layer  spacing  was  chosen  such  that  wall  function  spacing  was 
present  on  the  body,  while  a  solve-to-wall  mesh  was  used  on  the  fins  and  canards.  To  achieve 
this  end,  a  first  cell  spacing  for  the  prism  layer  ( dy )  was  specified  at  3.5  X  10-3  mm  on  the  fins 
and  canards,  0.25mm  on  the  body,  and  0.2  mm  on  the  nose.  MIME  allowed  for  a  smooth 
transition  between  the  spacings. 


Fig.  2  Extent  of  outer  mesh  boundary  with  close-up  of  projectile  showing  locations  of  density  boxes 

3.1.3  Flow  Field  and  Boundary  Conditions 

All  computations  were  completed  using  a  free-stream  temperature  and  pressure  of  288.15  K  and 
101325  Pa,  respectively.  During  the  course  of  the  study,  up  to  3  Mach  numbers  were 
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investigated:  Mach  0.50,  0.65,  and  0.80.  Up  to  14  angles  of  attack  between  —14°  and  14°, 
inclusive,  were  also  considered.  The  domain  was  initialized  using  free-stream  conditions 
everywhere.  For  some  cases,  the  area  immediately  surrounding  the  projectile  was  initialized  with 
zero  velocity  to  aid  convergence. 

The  entire  far- field  boundary  was  set  as  “characteristic  based”  inflow/outflow.  This  boundary 
condition  takes  the  specified  free-stream  conditions  and  solves  a  Riemann  problem  at  the 
boundary  using  the  supplied  data  as  a  virtual  state  outside  the  domain.  The  walls  of  the  projectile 
were  specified  as  adiabatic,  no  slip,  viscous  walls.  The  wall  function  integration  option  was 
chosen  for  the  projectile  surface  unless  the  first  prism  layer  was  such  to  ensure  proper  integration 
to  the  wall,  as  was  the  case  on  the  fins  and  canards. 

3.2  Parameter  Sensitivity  Technique 

The  goal  of  the  parameter  sensitivity  routine  was  to  assess  the  maneuver  characteristics  of 
various  candidates.  Airframe  input  data  (aerodynamics  and  mass  properties)  were  perturbed  to 
arrive  at  a  hypothetically  optimal  configuration.  Lift-to-drag  ratio  at  the  trim  angle-of-attack  was 
the  key  metric  for  optimization.  Canard  size,  canard  location,  canard  deflection,  fin  size,  fin 
location,  and  center-of-gravity  location  were  parameterized  from  the  nominal  values  to  derive 
theoretical  configurations. 

The  approach  for  the  parameter  sensitivity  routine  was  to  begin  by  assembling  the  baseline  mass 
properties  and  aerodynamics.  The  body,  fin,  and  canard  aerodynamics  were  considered 
separately  to  facilitate  the  optimization  analysis  and  to  allow  the  independent  parameters  of 
interest  to  be  varied.  These  independent  variables  included  canard  and  fin  size  (area),  canard 
deflection  angle,  and  canard  hinge  location.  Once  the  baseline  aerodynamic  coefficients  were 
specified,  the  canard  and  fin  size  could  be  altered  through  scaling  of  the  axial  and  normal  force 
coefficients.  Additionally,  the  canard  deflection  angles,  S,  were  accounted  for  by  adding  Sto  a 
and  determining  the  canard  aerodynamics  at  the  total  canard  angle.  Finally,  if  canard  location  or 
center  of  gravity  location  was  varied,  the  variation  was  accounted  for  by  properly  adjusting  the 
pitching  moment  coefficient. 

The  sensitivity  algorithm  uses  a  given  set  of  component  aerodynamic  coefficient  derivatives  as 
input  for  the  airframe.  Additionally,  a  range  of  values  for  each  parameter  is  specified  to 
determine  if  a  better  configuration  can  be  designed  within  the  limitations  of  the  given 
aerodynamic  coefficient  database.  The  only  constraint  on  the  optimization  is  that  a  given 
configuration  is  stable  (marginal  stability  is  allowed).  For  each  set  of  parameters,  the  pitching 
moment  of  the  individual  components  are  determined  and  then  summed  to  obtain  the  total 
pitching  moment  for  the  parameter  set  as  shown  in  Eq.  2. 
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(2) 


^  Cm(M,  a )  =  fFCFm(M,  a)  +  C«  (Af,  a)  +  /c[C£(Af,  5  +  a) 

+  (*C  -  *cre/)  Cw(M,  5  +  a)] 

In  this  equation,  fc  and  fF  are  the  canard  and  fin  scaling  factor,  respectively;  xc  is  the  canard 
hinge  location  forward  of  the  Xcg ;  xCref  is  the  canard  hinge  location  at  which  the  canard 
aerodynamic  coefficients  were  determined;  q  is  the  dynamic  pressure;  5  is  the  aerodynamic 
reference  area;  D  is  the  reference  diameter;  Cm  is  the  pitching  moment;  and  the  superscripts  F, 

B,  and  C  represent  the  fin,  body,  and  canard,  respectively.  The  trim  angle-of-attack  (atrim)  was 
found  by  finding  the  angle  for  which  Eq.  2  equals  zero  (i.e.,  Cm  changes  sign).  If  no  atrim  exists, 
the  parameter  set  is  considered  to  give  an  unstable  design  and  therefore  is  not  considered  further. 

The  total  axial  and  normal  force,  Cx  and  CN,  respectively,  at  this  atrim  was  determined  as  shown 
in  Eqs.  3  to  4  below. 

Cn(M,  atrim)  ,  dtrirn)  "h  (.M,  C^trirn)  "h  fc^N^^’8  "h  trim ) 

Cy(M,  d-trim)  /f^X  3"  (AC  (X-triTn)  3"  fc^X  (AC  8  +  ttfrim) 

Finally,  the  lift-to-drag  ratio  ( L/D )  was  calculated  from  Eq.  5. 

Y  j  ,  CCfrim)  COS((Z trim )  ^xC^C  ttfrim)  sin(atrim) 

,  Ctirim)  sia(£rtrim)  +  Cx (AC  OCtrim)  COS {d-trirn) 


(3) 

(4) 

(5) 


The  optimal  hypothetical  design  is  determined  by  finding  the  combination  of  parameters  that 
produce  the  maximum  L/D  at  its  dtrim.  This  sensitivity  analysis  is  relatively  general  and  can  be 
applied  to  various  aerodynamic  surfaces  with  different  aerodynamics  models. 

3.3  Flight  Mechanics  and  Control 

The  controlled  flight  behavior  of  the  different  airframe  configurations  was  investigated.  Models 
of  the  flight,  actuation,  and  measurements  were  used  in  flight  control  algorithms  and 
implemented  in  simulation  to  assess  performance.  A  linear  state  space  representation  (i  =  Ax  + 
Bu  +  F )  of  the  projectile  flight  with  a  first  order  model  for  the  actuator  is  presented.  The  state 
vector  is  composed  of  the  roll  angle  (0);  angular  rates  (p,  q,  r);  lateral  accelerations  (t>,  w);  and 
deflections  in  roll  ( 8p ),  pitch  (8q),  and  yaw  (8r)  as  shown  in  Eq.  6. 

X  =  [(f)  p  q  r  V  w  Sp  sq  8r]  (6) 
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The  controls  are  the  commanded  deflections  in  roll  (SCp),  pitch  (5Cq),  and  yaw  (SCr)  as  shown 
in  Eq.  7. 


U  —  \8c,p  ^C,q  SCx] 

The  system  dynamics  matrix,  with  mass  properties,  aerodynamic  characteristics,  and  actuator 
parameters,  is  provided  in  Eq.  8  below 


QSD  D 
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The  following  expression  for  controls  matrix  (Eq.  9)  contains  the  actuator  time  constant. 


The  aerodynamic  roll  torque  (Eq.  10)  is  the  only  term  in  the  F  matrix. 


(10) 


Feedback  of  roll  angle,  angular  rate  sensors,  and  accelerometers  was  utilized  in  this  guidance, 
navigation,  and  control  scheme,  which  yields  the  measurement  matrix  in  Eq.  11. 

^6x6  ^6X3 

^3X6  O3X3 

This  linear  model  was  implemented  in  simulation  and  then  used  with  optimal  control  theory  for 
preliminary  flight  control  design  of  the  high  maneuverability  airframe. 

The  nonlinear  flight  mechanics  of  these  airframes  was  also  studied.  The  6  degree-of-freedom 
model  is  provided.  The  translational  and  rotational  kinematic  equations  are  shown  in  Eqs.  12  and 
13,  respectively.  Shorthand  is  used  for  trigonometric  quantities  in  the  equations. 


ce  cip 

S(f>  Sq  C(p  $ip 

^•(fiSgC^,  T  S(p  Sxp 
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C(p/ce J  '-r J 

The  dynamic  equations  are  given  in  Eqs.  14  and  15.  The  forces  ( X ,  Y,  Z)  and  moments  (L,  M,  N ) 
acting  on  the  projectile  include  both  the  aerodynamic  and  gravitational  components. 

til  1  m  0  -r  q  1  it  ■ 

t)=  —  Y  —  r  0  —pv  (14) 

771 

-Wi  IZl  l—q  p  0  J  LwJ 

p]  fLl  [0  -r  q  I  rp- 

q  =  r1  M  - 1'1  r  0  -p  I  q  (15) 

.r  J  IN  J  -q  0  J  Lr- 

The  nonlinear  models  were  implemented  in  simulation  and  enabled  comprehensive  flight  control 
design.  The  flight  mechanics  and  control  for  the  high  maneuverability  airframe  are  provided  in 
greater  depth  in  Fresconi  et  al.16 

4.  Results  and  Discussion 

The  study  presented  herein  was  completed  in  2  phases.  The  first  phase  of  the  study  was  carried 
out  for  the  initial  body-fin-canard  configurations  that  were  designed  based  on  the  results  of  the 
component  build-up  in  Silton  and  Fresconi4  (i.e.,  aerodynamic  coefficients  of  the  body,  fins,  and 
canard  components  were  determined  separately  and  then  summed  for  the  total  aerodynamic 
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coefficients).  In  this  phase  of  the  study,  CFD  simulations  were  completed  to  determine  the  effect 
of  the  fin-canard  interactions  on  the  aerodynamic  coefficients  and,  therefore,  on  the  stability  of 
the  configuration,  as  well  as  on  munitions  in  general.  The  aerodynamic  results  of  the  first  phase 
were  used  as  the  aerodynamic  model  for  the  follow-on  parameter  sensitivity  study.  The  resulting 
parameter  set,  to  include  fin  and  canard  form  factors  and  canard  hinge  location,  was  used  as  the 
basis  for  an  exterior  redesign. 

In  the  second  phase,  CFD  simulations  were  completed  for  the  redesigned  configuration,  the 
results  of  which  were  used  to  further  refine  the  aerodynamic  model.  Model  refinement  included, 
but  was  not  limited  to,  investigating  the  use  of  higher  order  ballistic  fits  for  the  axial  force  and 
roll  torque  coefficients  for  each  of  the  components  (i.e.,  body,  fins,  individual  canards)  to 
determine  if  symmetry  about  a  =  0°  was  appropriate,  as  this  was  not  done  in  the  first  phase 
because  lower  order  polynomial  fits  (i.e.,  symmetry  was  not  enforced)  were  used.  The  refined 
aerodynamic  model  was  first  used  for  an  additional  parameter  sensitivity  study  to  ensure  that 
there  wasn’t  a  better  parameter  set  available.  Finally,  this  refined  aerodynamic  model  was  used 
in  a  flight  mechanics  and  control  analysis  to  verify  that  the  redesigned  configuration  would 
indeed  meet  the  requirements  of  the  demonstration  platform. 

4.1  Phase  1 

The  main  purpose  of  the  first  phase  of  the  study  was  to  account  for  the  effects  of  the  fin-canard 
interactions  on  the  aerodynamic  coefficients,  which  had  been  neglected  in  previous  studies.4 
Aerodynamic  computations  were  carried  out  based  on  designs  similar  to  those  in  Silton  and 
Fresconi4.  However,  the  fin-  and  canard-scaling  factors  found  to  be  optimal  in  Silton  and 
Fresconi4  could  not  be  implemented  due  to  physical  constraints  on  the  placement  and  size  of  the 
canards  (i.e.,  xc  =  16.6  mm  [0.2  cal.]  from  Xcg  was  not  achievable).  As  such,  a  slight  redesign 

was  performed  prior  to  beginning  this  phase. 

After  completing  a  trade  study,  the  decision  was  made  to  proceed  forward  with  2  configurations, 
Configuration  la  and  Configuration  lb.  The  center  of  gravity  location  (Xcg)  was  placed  at  200 
mm  (2.41  cal)  from  the  nose  to  be  consistent  with  the  previous  study.  Both  configurations  would 
have  4  rectangular  fins  with  a  22.0-mm  chord  and  188.4-mm  tip-to-tip  span,  corresponding  to 
approximately  75%  of  the  fin  area  used  in  Silton  and  Fresconi.4  A  double  bevel  geometry  was 
used  to  produce  roll  torque.  Both  configurations  also  had  4  NACA0015  canards,  which  were 
aligned  with  the  fins.  The  canard  chord  was  maintained  at  18.86  mm,  but  the  canard  span  and 
dither  axis  location  were  allowed  to  vary.  Configuration  la  had  a  tip-to-tip  canard  span  of  240.48 
mm  and  a  canard  axis  location  68.89  mm  (0.83  cal.)  forward  of  Xcg  (Fig.  3,  top).  Configuration 
lb  had  a  canard  axis  location  91.30  mm  (1.1  cal.)  forward  of  Xcg  and  a  192.72-mm  tip-to-tip  span 
(Fig.  3).  The  tip-to-tip  canard  span  for  the  canard  aerodynamic  data  used  for  the  parameter 
sensitivity  study  in  Silton  and  Fresconi4  was  approximately  191  mm,  so  that  the  canard  area  in 
that  study  was  similar  (-100%  canard  area,  fc  =  1.0)  to  the  canard  area  in  Configuration  lb. 
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Fig.  3  Initial  body-fin-canard  geometries  Configuration  la  (top)  and  Configuration  lb  (bottom) 

4.1.1  Static  Aerodynamic  Characterization 

The  static  aerodynamic  characterization  of  both  initial  configurations  was  completed  for  each  of 
the  3  Mach  numbers  (Mach  0.5,  0.65,  0.8)  at  zero  canard  deflection  angle.  Seven  angles-of- 
attack  (positive  angles  only)  were  investigated  for  Configuration  la  (21  simulations)  and  4 
angles-of-attack  (positive  angles  up  to  10°)  were  investigated  for  Configuration  lb  (12 
simulations).  For  each  simulation,  a  suitable  number  of  iterations  were  completed  (typically  2000 
to  3000)  to  ensure  convergence  in  both  residual  drop  as  well  as  in  aerodynamic  coefficients.  The 
aerodynamic  coefficients  were  obtained  by  averaging  their  values  over  the  last  200  iterations  to 
minimize  any  oscillations  (mostly  present  at  higher  angles  of  attack).  The  simulations  were 
completed  on  the  IBM  iDATAPlex  system,  Pershing,  at  the  ARL  DSRC.  Using  the  meshing 
technique  described  in  Section  3.1.2,  a  mesh  containing  approximately  33  million  cells  was 
created  for  both  Configurations  la  and  lb.  Each  33  million  cell  mesh  was  decomposed  using 
pMETIS  to  run  on  112  processors  (approximately  295,000  cells  per  processor),  which  resulted  in 
an  average  processing  time  of  approximately  13  s  per  iteration. 
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For  each  configuration,  the  axial  force  and  roll  torque  coefficients  ( Cx  and  Ch  respectively)  as 
well  as  normal  force  and  pitching  moment  coefficients  ( CN  and  Cm,  respectively)  were 
determined  at  each  Mach  number  and  each  angle  of  attack.  The  forces  and  moments  were 
examined  for  the  complete  configuration  as  well  as  for  the  contributions  of  the  individual 
components  (i.e.,  body,  total  fins,  and  each  individual  canard).  The  coefficient  data  for  each  of 
the  components  was  transformed  from  the  CFD  coordinate  system  to  the  projectile  body 
coordinate  system.  The  coefficient  data  for  the  canards  was  further  converted  to  the  local  blade 
coordinate  system  (the  y-axis  aligned  with  the  blade,  the  z-axis  perpendicular  to  the  blade).  The 
relative  orientation  of  the  canards  is  shown  in  Fig.  4;  the  roll  orientation  for  Canard  4  is  45°  with 
the  angle  increasing  in  the  clockwise  direction  (i.e.,  decreasing  canard  number). 


Fig.  4  Canard  location  looking  from  rear  of  projectile 

A  detailed  explanation  of  these  coordinate  systems  and  coordinate  transformations  can  be  found 
in  Appendix  A.  Polynomial  and  ballistic  fits  of  the  coefficients  (Eqs.  16  tol9)  were  then 


determined  at  each  Mach  number  for  each  component. 

Cx  =  CXo  +  CXal8  +  CXa28 2  +  CXa3S 3  +  CXa4S4  (16) 

CN  =  CNaS  +  CNa383  (17) 

Cm  =  Cma8  +  Cma3S 3  (18) 

Ci  =  Clo  +  Clal8  +  Cla282  +  Cla3S3,  (19) 

where 

8  =  sin  a.  (20) 


Figures  5  and  6  show  the  ability  of  the  fits  to  capture  the  data  trends  at  Mach  0.5  for  the  body 
and  fin  components  and  the  individual  canard  components,  respectively,  for  Configuration  la. 
The  fits  were  equally  as  good  for  the  other  Mach  numbers  as  well  as  for  Configuration  lb.  A 
complete  set  of  plots  showing  the  coefficient  fits  can  be  found  in  Appendix  B.  It  is  important  to 
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remember  that  only  positive  body  angles  of  attack  were  considered  in  the  first  phase. 
Polynomial,  rather  than  even  ballistic  fits,  were  implemented  for  Cx  and  Ct  to  ensure  the  best  fit 
to  the  data.  However,  the  use  of  polynomial  fits  does  not  guarantee  symmetry  about  a  =  0°, 
which  may  be  more  realistic. 


Fig.  5  Body  and  Fin  data  and  fits  in  projectile  body  coordinate  system  for  Configuration  la,  Mach  0.5 


Ct  is  fit  for  both  the  total  fins  and  individual  canards  (body  contribution  is  zero).  However,  the 
total  canard  contribution  in  the  CFD  simulations,  as  the  canards  are  not  deflected,  is  zero.  The 
resulting  fits  for  each  component  are  converted  to  ballistic  coordinates.  Tables  2  and  3  list  the 
values  of  coefficient  derivatives  in  ballistic  coordinates  for  both  configurations  for  the  body  and 
fins  components,  respectively.  For  the  fin  component,  only  a  second  order  ballistic  fit  was 
required  to  accurately  fit  the  axial  force  coefficient,  Cx;  only  CXq  and  Cx  are  listed  in  Table  3  as 

the  remainder  of  the  coefficients  have  a  value  of  zero. 

Each  of  the  canard  aerodynamic  fits  was  determined  independently  in  the  local  canard  coordinate 
system.  This  means  Canards  2  and  3  have  a  negative  local  blade  angle  and  Canards  1  and  4  have 
a  positive  local  blade  angle.  As  with  the  body  and  fin  data,  polynomial  fits,  rather  than  ballistic 
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fits,  were  implemented  for  Cx  and  Cl  to  ensure  the  best  fit  to  the  data  rather  than  forcing 
symmetry  about  a  =  0°.  Tables  4  to7  list  the  values  of  coefficient  derivatives  in  ballistic 
coordinates  (local  canard  coordinate  system)  for  both  configurations  for  each  canard, 
respectively.  There  is  no  value  listed  for  CLq  as  it  is  always  zero  due  to  the  canards  not  being 

deflected.  These  coefficient  derivatives  are  used  in  Eqs.  16  tol9  to  determine  the  value  of  the 
coefficient  attributable  to  each  component  at  each  Mach  number  for  a  given  angle  of  attack. 


o* 

<L> 

-o 

i? 

CD 

15 

o 

o 


Local  Blade  Angle-of-Attack  (deg) 


-10  -5  0  5  10 

Local  Blade  Angle-of-Attack  (deg) 


oN 

<D 

TD 

CD 

15 

o 

o 


-0.25 


•  Canard  #1 

- Canard  #1  (Fit) 

- 

■  Canard  #2 

- Canard  #2  (Fit) 

- 

♦  Canard  #3 

- Canard  #3  (Fit) 

~ 

A  Canard  #4 

Canard  #4  (Fit) 

0.2 


0.1 


—  n 


-0.1 


-10  -5  0  5  10 

Local  Blade  Angle-of-Attack  (deg) 


-0 


•3, 


•  Canard  #1 

- Canard  #1  (Fit) 

■  Canard  #2 

- Canard  #2  (Fit) 

♦  Canard  #3 

- Canard  #3  (Fit) 

A  Canard  #4 

. Canard  #4  (Fit) 


-10  -5  0  5  10 

Local  Blade  Angle-of-Attack  (deg) 


15 


Fig.  6  CFD  data  and  fits  in  for  each  canard  in  local  canard  coordinate  system  for  Configuration  la,  Mach  0.5 

The  component  coefficients  are  used  as  input  for  the  flight  dynamics  model.  During 
implementation,  the  canard  contribution  must  be  transformed  from  the  local  blade  coordinate 
system  back  to  the  projectile  body  coordinate  system.  Figure  7  shows  the  results  of  this 
transformation  for  Configuration  la  at  Mach  0.5.  Because  the  canards  are  in  an  orientation  that  is 
45°  off  plane,  the  canard  angle  of  attack  in  the  local  blade  coordinate  system  is  smaller  than  that 
in  the  projectile  body  coordinate  system.  Additionally,  only  a  portion  of  the  force  along  the 
canard  z-axis  is  the  Z-force  direction  of  the  projectile  body.  Therefore,  the  slopes  of  the 
coefficients  are  decreased  during  the  transformation  from  the  local  blade  coordinate  system  back 
to  the  body  coordinate  system. 
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Table  2  Body  component  of  aerodynamic  coefficient  fits  for  Configurations  la  and  lb 


Config. 

Mach 

C*o 

CXal 

CXa2 

C*«3 

C  ]y 

7V  a 

CN«3 

r 

r 

^rna3 

la 

0.50 

0.1506 

0.1282 

-0.7824 

3.0321 

3.8723 

-13.2237 

4.1850 

-4.1115 

0.65 

0.1482 

-0.0224 

1.3919 

-3.1198 

3.6772 

-13.6664 

4.2459 

-5.8507 

0.80 

0.2128 

0.0124 

2.4050 

-4.9798 

2.7188 

3.0697 

3.9527 

-2.6088 

lb 

0.50 

0.1529 

0.0000 

0.3501 

0.0000 

3.6358 

-10.9442 

3.77247 

12.3314 

0.65 

0.1482 

0.0000 

0.7685 

0.0000 

3.8299 

-24.9079 

4.1625 

-2.9340 

0.80 

0.2140 

0.0000 

1.5736 

0.0000 

2.9240 

-7.5938 

3.9940 

-3.7350 
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Fig.  7  Comparison  of  canard  coefficients  in  local  blade  coordinate  and  projectile  body  coordinates  for 
Configuration  la,  Mach  0.5 
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Table  3  Total  fin  component  of  aerodynamic  coefficient  fits  for  Configurations  la  and  lb 


Config. 

Mach 

CXo 

CX«2 

C  ]y 

CNa  3 

r 

L’rna 

r 

C,0 

cial 

Cia2 

Cia, 

la 

0.50 

0.05413 

0.31131 

1.75757 

-6.0255 

-3.8427 

12.8555 

0.1143 

0.2975 

-5.6023 

12.2867 

0.65 

0.05536 

0.39140 

1.72327 

-5.3884 

-3.7730 

11.4362 

0.1196 

0.2860 

-6.6319 

16.1398 

0.80 

0.05662 

0.55466 

1.74705 

-5.6446 

-3.8301 

11.9842 

0.1200 

-0.1006 

-3.6367 

10.9799 

lb 

0.50 

0.05424 

0.27124 

1.9314 

-10.5717 

-4.2027 

21.9989 

0.1148 

0.1578 

-2.9050 

0.37678 

0.65 

0.05485 

0.43284 

1.85519 

-8.7736 

-4.0443 

18.1562 

0.1198 

0.2025 

-5.1952 

9.9199 

0.80 

0.05484 

0.75421 

1.82664 

-7.9024 

-4.0004 

16.6889 

0.1202 

0.06957 

-6.9581 

24.7248 

Table  4 

Aerodynamic  coefficient  fits  Canard  lfor  Configurations  la  and  lb 

Config. 

Mach 

cx0 

CX„1 

CXa2 

CXa3 

C*«4 

C n 

CJV«3 

Cma 

Cma  3 

Clal 

Cla  2 

Cla  3 

la 

0.50 

0.00614 

0.0776 

-3.879 

14.984 

0.0 

1.953 

-22.118 

1.638 

-18.555 

-1.474 

— 4.2 14 

32.995 

0.65 

0.00616 

0.0712 

-5.558 

52.944 

-137.508 

1.845 

-32.249 

1.561 

-27.530 

-2.165 

6.602 

4.458 

0.80 

0.0121 

-0.00589 

0.760 

-2.102 

-4.745 

0.817 

-2.991 

0.708 

-3.413 

-0.832 

1.273 

-2.724 

lb 

0.50 

0.00417 

0.0 

-0.941 

0.0 

0.0 

1.134 

-5.254 

1.291 

-5.842 

-0.880 

-0.106 

4.247 

0.65 

0.00410 

0.0182 

-2.178 

13.730 

0.0 

1.339 

-32.392 

1.530 

-36.963 

-0.909 

-3.220 

41.574 

0.80 

0.00718 

0.00287 

0.177 

0.928 

0.0 

0.709 

-14.676 

0.826 

-17.460 

-0.630 

2.331 

-3.268 

Table  5  Aerodynamic  coefficient  fits  Canard  2  for  Configurations  la  and  lb 


Config. 

Mach 

Cx0 

Cxai 

CXa2 

Cxa3 

CXai 

CNa 

CNa  3 

Cma 

r 

^™a3 

Cla  1 

Cla  2 

Cla3 

0.50 

0.00622 

0.0650 

-3.622 

-13.832 

0.0 

1.934 

-21.547 

1.621 

-18.058 

-1.503 

3.616 

30.426 

la 

0.65 

0.00617 

-0.0701 

-5.549 

-53.139 

-138.790 

1.837 

-32.469 

1.555 

-27.689 

-2.177 

-6.852 

3.775 

0.80 

0.0121 

0.00526 

0.736 

1.909 

-5.338 

0.777 

-2.302 

0.678 

-2.892 

-0.810 

-1.425 

-3.775 

0.50 

0.00416 

0.0 

-0.947 

0.0 

0.0 

1.134 

-5.152 

1.291 

-5.724 

-0.880 

0.113 

4.215 

lb 

0.65 

0.00410 

-0.0173 

-2.153 

-13.505 

0.0 

1.339 

-32.640 

1.530 

-37.196 

-0.909 

3.223 

41.834 

0.80 

0.00718 

-0.00316 

0.103 

-1.356 

0.0 

0.693 

-14.195 

0.809 

-16.923 

-0.611 

-2.164 

-  -2.623 

Table  6  Aerodynamic  coefficient  fits  Canard  3  for  Configurations  la  and  lb 


Config. 

Mach 

Cx0 

CXal 

CXa2 

Cxa3 

Cxai 

C  ft 

CN„3 

Cma 

Cma  3 

Ca 

0.50 

0.00630 

-0.0468 

-3.217 

-12.367 

0.0 

1.837 

-18.004 

1.540 

-15.278 

-1.578 

1.318 

18.452 

la 

0.65 

0.00616 

-0.0725 

-5.716 

-55.545 

-148.299 

1.773 

-27.253 

1.498 

-23.442 

-2.127 

-6.714 

-1.650 

0.80 

0.0121 

-0.00467 

1.043 

7.837 

16.518 

0.814 

1.451 

0.702 

0.154 

-0.667 

2.104 

5.684 

0.50 

0.00410 

0.0 

-0.882 

0.0 

0.0 

1.128 

-6.360 

1.283 

-7.208 

-0.879 

0.0574 

4.443 

lb 

0.65 

0.00410 

-0.0173 

-2.175 

-13.954 

0.0 

1.322 

-31.225 

1.510 

-35.773 

-0.916 

2.717 

36.639 

0.80 

0.00718 

-0.00914 

0.657 

3.425 

0.0 

0.573 

-1.833 

0.672 

-3.218 

-0.575 

-2.789 

-15.596 

Table  7  Aerodynamic  coefficient  fits  Canard  4  for  Configurations  la  and  lb 


Config. 

Mach 

Cx0 

cx , 

Aa  1 

Cxa2 

Cxa3 

CX  A. 

C n 

CNU3 

r 

^ma 

r 

Clal 

Cla2 

Cla  3 

0.50 

0.00628 

0.0514 

-3.307 

12.763 

0.0 

1.841 

-18.032 

1.544 

-15.320 

-1.574 

-1.418 

18.812 

la 

0.65 

0.00617 

0.0708 

-5.640 

54.752 

-145.926 

1.772 

-26.793 

1.497 

-23.067 

-2.128 

6.760 

-2.299 

0.80 

0.0121 

0.00283 

1.134 

-8.612 

18.497 

0.841 

1.245 

0.723 

0.0124 

-0.701 

-1.891 

4.894 

0.50 

0.00410 

0.0 

-0.874 

0.0 

0.0 

1.129 

-6.577 

1.284 

-7.459 

-0.876 

-0.135 

5.013 

lb 

0.65 

0.00410 

0.183 

-2.204 

14.192 

0.0 

1.325 

-31.317 

1.514 

-35.906 

-0.914 

-2.797 

37.132 

0.80 

0.00718 

0.00968 

0.653 

-3.381 

0.0 

0.597 

-2.485 

0.698 

-3.927 

-0.581 

2.521 

-13.786 

Once  the  component  contributions  have  been  determined,  the  total  aerodynamic  forces  and 
moments  can  be  determined  by  adding  these  contributions  together  (i.e.,  component  build  up). 
For  this  study,  the  component  contributions  include  the  interference  effects  between  the 
components;  they  are  not  considered  separately.  As  can  be  seen  from  Fig.  8  (Configuration  la, 
Mach  0.65)  and  Fig.  9  (Configuration  lb,  Mach  0.65),  when  the  component  fits  are  added 
together,  there  is  agreement  with  the  total  coefficients  obtained  from  CFD,  which  validates  that 
the  correct  component  fits  are  being  used.  (All  of  component  fit  build  ups  for  Configurations  la 
and  lb  can  be  found  in  Appendix  B).  Note,  CA  is  used  for  axial  force  coefficient  in  these  figures, 
rather  than  Cx,  as  the  aerodynamic  forces  are  being  represented  as  the  ballistic  coefficients. 
However,  it  is  quite  apparent  that  neither  of  these  configurations  is  statically  stable  (i.e.,  Cm  >  0) 
as  was  predicted  during  the  trade  study  completed  while  designing  the  projectile.  The  trade  study 
used  the  coefficients  obtained  from  body-fin  simulations  and  the  previously  obtained  canard  data 
presented  in  Silton  and  Fresconi.4  The  question,  of  course,  is  why  the  trade  study  predicted 
different  results  than  were  determined  by  the  CFD  simulations. 


Fig.  8  Component  fit  build  up  comparison  to  total  CFD  coefficients  for  Configuration  la,  Mach  0.65 
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Fig.  9  Component  fit  build  up  comparison  to  total  CFD  coefficients  for  Configuration  lb,  Mach  0.65 

The  lack  of  airframe  stability  can  be  answered  by  comparing  the  component  aerodynamic 
contributions  found  during  the  trade  study  to  those  component  aerodynamic  contributions 
determined  directly  from  the  CFD  (Fig.  9);  specifically  the  body  and  fin  components.  The  results 
of  the  trade  study,  which  used  the  appropriately  scaled  body,  fin,  and  canard  aerodynamic 
coefficients  (75%  fin  effectiveness,  100%  canard  effectiveness)  from  Silton  and  Fresconi4  to 
correspond  to  Configuration  lb,  are  shown  in  Fig.  10.  No  canard-fin  interaction  effects  are 
included  in  the  trade  study.  The  differences  between  the  aerodynamic  contributions  of  the  body 
and  fins  are  quite  apparent.  At  a  =  2°,  there  was  less  fin  force  (20%)  and  moment  (45%) 
predicted  by  the  CFD  for  the  body-fin-canard  model  than  the  body-fin  model  (Fig.  11),  which 
can  likely  be  attributed  to  the  interaction  effects  of  the  canard  tip  vortices  with  the  fin  blades. 
Additionally,  there  are  increases  in  body  normal  force  (20%)  and  pitching  moment  (60%)  (see 
Fig.  12),  which  are  also  likely  due  to  the  presence  of  the  canards.  These  results  show  just  how 
important  it  is  that  interaction  effects  be  accounted  for  when  creating  the  aerodynamic  model.  It 
is  obvious  that  the  stability  of  the  projectile  may  not  be  as  expected  if  the  canard-fin  or  the 
canard-body  interaction  effect  is  neglected.  Finally,  the  canards  themselves  show  a  maximum  of 
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a  25%  increase  in  normal  force  and  pitching  moments  at  angles  of  attack  below  7°  (Fig.  13), 
which  caused  the  projectile  to  be  less  stable  than  expected. 

This  large  of  a  variation  in  the  canard  contribution  over  the  experimental  values  used  in  Silton 
and  Fresconi4  was  a  bit  surprising;  the  same  interaction  effects  causing  the  increase  in  body 
normal  force  and  pitching  moment  may  be  causing  the  increased  canard  effectiveness.  It  is  also 
possible  that  the  different  body  shape  (ogival  rather  than  cylindrical)  on  which  the  canards  were 
mounted  in  the  experiment  is  causing  the  differences  in  canard  aerodynamic  forces.  The 
differences  become  more  significant  above  a  =  7°,  as  the  slope  of  the  canard  component  from 
the  CFD  appears  to  be  decreasing,  which  indicates  that  aerodynamic  stall  may  be  beginning  to 
occur,  while  the  experimental  data  used  in  the  component  build-up  indicates  that  aerodynamic 
stall  does  not  occur  until  nearly  a  =  12°.  Wind  tunnel  tests  will  be  conducted  in  the  future  to 
investigate  the  effect  of  the  body  shape  and  the  onset  of  aerodynamic  stall. 
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Fig.  10  Component  build  up  of  normal  force  and  pitching  moment  using  body-fin  CFD  and  canard  fit  data 
from  Silton  and  Fresconi  with  appropriate  form  factors  to  approximate  Configuration  lb,  Mach 
0.65  (fF=  0.75,  fc=  1.0) 


21 


Fig.  1 1  Comparison  of  fin  component  as  determined  from  component  buildup  with  form  factors  and  CFD 
simulation  for  Configuration  lb,  Mach  0.65. 


Fig.  12  Comparison  of  body  component  as  determined  from  component  buildup  with  form  factors  and  CFD 
simulation  for  Configuration  lb,  Mach  0.65 
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Fig.  13  Comparison  of  total  canard  component  as  determined  from  component  buildup  with  form  factors  and  CFD 
simulation  for  Configuration  lb,  Mach  0.65 

4.1.2  Redesign  for  Static  Stability  and  Parameter  Sensitivity  Study 

The  aerodynamic  data  for  both  Configurations  la  and  lb  were  used  as  input  to  the  numerical 
implementation  of  Eqs.  2  to5  to  redesign  the  airframe  to  ensure  static  stability  while  maximizing 
lift-to-drag  ratio.  The  goal  was  to  apply  the  results  of  a  high  fidelity  aerodynamic 
characterization,  which  included  the  flow  interaction  effects,  to  obtain  attractive  airframe 
specifications.  A  sensitivity  study  was  performed  to  determine  the  responsiveness  of  the  designs 
to  variations  in  fin  size  and  canard  size  and  location.  Canard  deflection  angle  was  also  varied  to 
determine  the  trim  angle  that  would  maximize  the  lift-to-drag  ratio.  Because  the  fin 
aerodynamics  included  interaction  effects  due  to  the  canards,  scaling  the  fin  area  (through  fF) 
scaled  the  interaction  effects  as  well. 

Aerodynamic  data  were  only  available  for  the  non-deflected  canard  case.  However,  data  for  the 
deflected  canard  case  were  necessary  to  use  the  equations.  Because  the  canard  coefficient  fits 
were  determined  in  the  local  blade  coordinate  system,  the  deflection  angle  is  added  to  the  local 
angle  of  attack;  this  total  angle  of  attack  is  then  used  to  determine  the  aerodynamic  coefficient  of 
the  canard.  Figure  14  shows  the  results  of  this  methodology  to  determine  the  total  coefficients 
assuming  a  deflection  angle  of  4°  for  a  pitch  maneuver  (i.e.,  all  deflected  “up”)  using 
Configuration  la.  While  this  method  is  not  exact,  it  is  a  reasonable  assumption  at  this  time. 
Future  studies  will  investigate  its  validity.  Additionally,  canard  stall  was  accounted  for  by 
holding  the  canard  normal  force  and  pitching  moment  constant  for  total  angles  of  attack  outside 
the  range  for  which  CFD  was  completed,  which  is  the  reason  the  values  of  CN  and  Cm  remain 
constant  above  a  body  angle  of  attack  of  approximately  8°. 
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Fig.  14  Component  fit  build  up  of  normal  force  and  pitching  moment  with  4° 
canard  deflection  for  Configuration  la,  Mach  0.65 

Configuration  la  is  statically  unstable,  which  necessitates  larger  fins  and  smaller  canards.  As 
such,  the  parameter  sensitivity  routine  was  executed  using  the  aerodynamic  coefficients  obtained 
from  CFD  simulations  for  Configuration  la  scaled  by  values  of  fF  >  1  and  fc  <  1,  while 
varying  the  canard  deflection  angle.  The  overall  effects  on  lift  and  drag  with  nonlinearities  in 
angle-of- attack  produces  a  maximum  lift-to-drag  around  fF  =  1.95  for  fc  =  0.75  and  5  =  4°. 
Figure  15  shows  the  results  of  the  optimization  routine  at  these  conditions.  Each  plot  features  the 
trends  in  all  3  metrics,  although  the  desire  is  to  maximize  on  the  lift-to-drag  ratio  at  its  trim  angle 
of  attack,  as  a  function  of  a  specific  variable  with  all  other  variables  held  constant  at  the  optimal 
value.  For  example,  the  top  left  plot  shows  the  variation  in  trim  angle,  lift-to-drag,  and  normal 
force  as  a  function  of  fin  size  with  canard  area  scaled  by  0.75  and  the  canard  deflection  angle  set 
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at  4°.  The  maximum  lift-to-drag  ratio  of  3  occurs  for  a  fin  scaling  factor  of  1.95  although  the 
variation  is  not  large  for  the  scaling  factors  shown.  Reducing  the  size  of  the  canards  below  0.75 
has  a  much  greater  effect  on  reducing  the  lift-to-drag  ratio  (top  right  plot  of  Fig.  15).  Canard 
scaling  factors  above  0.75  are  not  shown  as  the  projectile  wasn’t  statically  stable.  The  measures 
of  performance  are  relatively  flat  with  canard  deflection  angle  due  to  canard  flow  separation 
effects. 


Fig.  15  Parameter  sensitivity  study  using  Configuration  la  aerodynamic  data  with//  =  1.95,  fc  =  0.75,  and  5  =  4° 

Figure  16  shows  some  optimization  results  using  the  aerodynamic  coefficients  obtained  for 
Configuration  lb  scaled  by  fF  =  1.7  and  fc  =  0.75  for  8  =  4°.  Overall,  the  trends  are  consistent 
with  the  analysis  for  Configuration  la  in  Fig.  15.  Configuration  lb  has  smaller  canards,  which 
means  that,  to  get  similar  metrics  as  shown  in  Fig.  15  for  Configuration  lb,  the  fin  size  does  not 
need  to  be  increased  as  much.  The  results  compiled  in  Figs.  15  and  16  suggest  that  a  larger  total 
fin  area  (i.e.,  more  and/or  larger  fins)  is  necessary.  The  exact  amount  will  be  determined  after 
completing  a  comprehensive  redesign  to  obtain  an  accurate  Xcg  location. 
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Fig.  16  Parameter  sensitivity  study  using  aerodynamic  data  from  Configuration  lb  with//?  =  1.7,  fc  =  0.75,  and 
6  =  4 0 

4.2  Phase  2 

The  main  purpose  of  the  second  phase  of  this  study  was  to  validate  that  the  optimized  airframe 
design  as  determined  in  the  first  phase  would  meet  the  requirements  to  be  used  as  the 
demonstration  platform  for  the  LCHAW  mission  program.  The  outcome  of  the  first  phase  was  an 
optimized  airframe  design  that  accounted  for  the  interaction  effects  between  the  body,  fins,  and 
canards  of  the  airframe.  Prior  to  validating  the  aerodynamic  performance  of  this  design  with 
additional  CFD  predictions,  efforts  were  undertaken  to  further  refine  the  design,  such  as  better 
determination  of  mass  properties,  especially  for  the  Xcg  location.  Figure  17  shows  the  airframe 
components  and  their  relative  assumed  placement,  which  were  included  in  the  solid  model  of  the 
projectile  to  calculate  mass  properties.  A  number  of  other  assumptions  in  designing  the  updated 
airframe  configuration  were  made,  including  that  the  volume  and  mass  of  the  potted  electronics 
and  batteries  be  fixed  and  that  they  be  positioned  directly  behind  the  control  actuation  system 
(CAS)  section,  which  was  located  directly  behind  the  camera.  It  was  assumed  that  the  canards 
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within  the  CAS  were  stowed  folded  forward.  Finally,  it  was  assumed  that  the  fins  were  stowed 
folded  forward  and  the  fin  hinge  point  was  fixed.  As  such,  the  deployed  fin  location  was  also 
fixed. 


8  Fins 


Batteries 


Potted 

Electronics 


4  Canards 


Camera 


Fig.  17  Projectile  components  used  to  calculate  mass  properties 


As  with  any  design,  there  were  a  number  of  constraints  that  had  to  be  met.  The  most  limiting 
constraint  was  that  the  overall  projectile  length  be  fixed  at  427.3  mm.  As  the  canards  are 
assumed  to  fold  forward  within  the  CAS,  the  forward  most  canard  pivot  location  needs  to  be 
constrained  by  the  canard  length  (i.e.,  shorter  canards  can  be  positioned  closer  to  the  nose  than 
long  canards).  Then,  as  fin  length  is  adjusted  for  stability,  a  corresponding  change  in  the  length 
of  the  fin  hub/boattail  section  needs  to  occur.  Therefore,  the  rearmost  limit  of  the  guidance 
electronics  section  and  the  CAS  section  locations  were  adjusted  by  the  same  amount.  This 
adjustment  constrains  the  rearmost  canard  location.  However,  this  constraint  was  not  a  limiting 
factor  for  any  of  the  configurations  that  were  evaluated  because  the  guidance  electronics  never 
impinged  on  the  fin  hub.  These  constraints  are  illustrated  in  Fig.  18. 

Using  these  additional  geometric  constraints,  the  parameter  sensitivity  study  was  completed  with 
the  aerodynamic  coefficients  of  Configuration  la  and  updated  Xcg.  By  parameterizing  fin  area, 
canard  area,  and  canard  axis  location,  Configuration  2  was  designed. 


27 


Large  Fins, 
Longer  Fin  Hub 


Fig.  18  Illustration  of  geometry  constraints 

Configuration  2  (Fig.  19)  has  8  rectangular  fins  with  a  199.87-mm  tip-to-tip  span,  a  19.05-mm 
chord  (200%  of  Configuration  la  fin  area),  and  a  1.27-mm  thickness  with  the  trailing  edge 
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located  5.10  mm  from  the  base  of  the  projectile.  For  ease  of  construction  in  the  future,  it  was 
decided  to  use  a  bent  fin  tip  to  create  roll  torque  rather  than  using  a  double  bevel  geometry  that 
had  been  used  in  Configurations  la  and  lb.  The  top  8.79  mm  of  each  fin  of  Configuration  2  was 
bent  to  an  angle  of  15°.  Four  NACA0015  airfoils  in  an  “X”  orientation  (i.e.,  45°  off  plane)  were 
used  for  the  canards.  The  canard  chord  was  18.86  mm  and  the  tip-to-tip  span  was  192.72  mm 
(70%  of  Configuration  la  canard  area).  The  canard  axis  (at  the  quarter-chord  of  the  canard)  was 
located  73.787-mm  (0.889  cal.)  forward  of  the  Xcg.  For  the  fin-canard  offset,  2  fins  were  equally 
spaced  between  each  set  of  canards  for  an  alignment  offset  of  22.5°.  Finally,  to  simulate  a 
possible  deployment  mechanism  for  the  fins  in  the  future,  a  cavity  was  present  in  the  base  of  the 
projectile  along  with  a  pin.  The  mass  properties  were  calculated  to  be  Xcg  =  218  mm  from  the  nose 
of  the  projectile,  Ix  =  0.0034  kg-m2,  and  Iy  =  0.0388  kg-m2,  and  m=  2.65  kg. 

4.2.1  Aerodynamic  Characterization 

The  static  aerodynamic  characterization  of  Configuration  2  was  completed  for  only  1  Mach 
number,  Mach  0.65,  as  the  change  in  mass  properties  indicated  that  the  velocity  would  not  vary 
significantly  during  its  flight.  Fourteen  angles  of  attack,  —14°  <  a  <  14°,  were  investigated.  For 
each  simulation,  sufficient  iterations  were  completed  (typically  2000  to3000)  to  ensure 
convergence  in  both  residual  drop  and  aerodynamic  coefficients.  The  aerodynamic  coefficients 
were  obtained  by  averaging  their  values  over  the  last  200  iterations  to  minimize  any  oscillations 
(mostly  present  at  higher  angles  of  attack).  The  simulations  were  again  completed  on  the  IBM 
iDATAPlex  system,  Pershing,  at  the  ARL  DSRC.  For  this  configuration,  a  mesh  containing 
approximately  40.6  million  cells  was  created.  The  40.6  million  cell  mesh  was  decomposed  using 
pMETIS  to  run  on  272  processors  (-149,000  cells/processor),  which  resulted  in  an  averaging 
processing  time  of  approximately  9  s  per  iteration. 

Only  ballistic  fits  were  required  for  the  coefficient  fits  for  the  body  and  fin  components  with  this 
configuration  as  inclusion  of  negative  angles  of  attack  here  prove  that  Cx  and  CL  are  in  fact 
symmetric  across  a  =  0°.  As  such,  the  fit  for  Cx  reduces  from  Eqs.  16  to  21  and  the  fit  for  Ct 
reduces  from  Eqs.  19  to  22.  The  component  contributions  to  the  aerodynamic  coefficients  for  the 
body  and  fins  (total)  are  listed  in  Table  8.  The  side  force  and  moment  contributions  are  currently 
considered  to  be  negligible.  However,  this  assumption  is  likely  not  to  be  the  case  if  and  when  the 
canards  are  not  all  deflected  equally.  Thus,  the  body  and  total  fin  side  force  and  moment 
contributions  will  have  to  be  revisited  as  the  design  matures.  Figure  20  shows  how  well  these  fits 
are  able  to  predict  the  aerodynamic  coefficients  as  determined  by  CFD.  The  fits  for  Cx  and  Ct 
predict  the  overall  shape  fairly  well  though  the  values  near  a  =  0°  are  not  as  well  predicted  as 
they  might  be.  At  this  point  in  the  design  process,  the  discrepancy  in  the  predictions  at  a  =  0°  is 
not  a  concern.  It  will  be  investigated  further  in  the  future. 

Cx  =  CX o  +  cxa2S2  +  CXa4S4  (21) 

ci  =  ci0  +  cia2s 2  +  cia4S4  (22) 
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Table  8  Static  aerodynamic  coefficient 
fits  for  Configuraion  2  at  Mach 
0.65,  body  and  fin  components 


C*o 

Body 

Fins 

0.1538 

0.1474 

CXa2 

-0.1943 

0.7101 

CXa* 

6.5204 

-4.8400 

CNa 

3.7278 

4.4948 

CNca 

-11.3042 

-25.6679 

Cma 

3.9785 

-10.4992 

r 

^™CC3 

-7.2300 

59.3082 

CU> 

0.08242 

Cla2 

0.6111 

Cla4 

-9.5560 
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Fig.  20  Body  and  fin  component  fits  for  Configuration  2,  Mach  0.65 

Figure  21  shows  the  fits  for  4  of  the  6  static  aerodynamic  coefficients  for  the  canards  in  the  local 
blade  coordinate  systems;  the  side  force  and  side  moment  are  assumed  to  be  0  in  the  local  blade 
coordinate  system.  Again,  because  the  canard  fits  are  in  the  local  blade  coordinate  system,  the 
side  force  and  moment  are  assumed  to  be  negligible.  The  numbering  scheme  for  the  canards 
remains  the  same  as  for  the  previous  configuration  (Fig.  4).  For  all  the  canard  aerodynamic 
coefficients,  it  is  quite  apparent  that  the  aerodynamic  behavior  changed  at  approximately 
a  =  ±5.5°;  Cz  and  Cm  were  approximately  constant  when  |a|  >  5.5°.  Thus,  it  is  likely  that 
aerodynamic  stall  of  the  canards  is  occurring  near  this  a.  As  such,  a  trimodal  fit  ( a  >  —5.5°,  a  > 
5.5°,  —5.5°  <  a  <  5.5°)  was  developed  for  the  canard  aerodynamic  coefficients  to  best  represent 
the  data  trends.  To  fit  the  data  around  a  =  0°,  only  ballistic  fits  were  required  for  the  individual 
canard  components.  As  such,  Eq.  21  is  again  used  for  Cx ■  However,  Ct  for  the  canards  is  an  odd 
fit  (the  constant  term  is  still  zero).  This  means  that  Eq.  19  reduces  to  Eq.  23  rather  than  Eq.  22. 

At  the  larger  absolute  angles  of  attack,  after  canard  stall  occurs,  the  coefficient  fits,  for  all 
aerodynamic  coefficients  of  interest,  change  drastically;  Cx  is  an  odd  fit  (Eq.  24),  while  Cn  (-Cz), 
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Cm,  and  Ci  are  linear  fits  (Eqs.  25  to  27,  respectively). To  ensure  that  the  transition  between  the 
segments  of  the  trimodal  fits  is  smooth,  the  value  of  the  constant  ( CX() ,  CNa0,  Cma0,  Cl(), 

respectively)  in  the  fits  beyond  stall  is  taken  as  the  value  of  the  coefficient  at  the  stall  angle,  as 
appropriate. 


C,  =  C,  S  +  C,  83 

1  la l  La3 

(23) 

=  CXo  +  CXal8  +  CXa3S 3 

(24) 

CN  =  CNa0  +  cNals 

(25) 

Cm  ~  Cma0  Cmal8 

(26) 

to 

s 

+ 

o 

O' 

II 

cr 

(27) 

The  component  contributions  of  each  of  the  4  individual  canards  for  each  of  the  a  ranges  are 
listed  in  Tables  9  to  12.  After  completing  the  fits,  a  more  exact  value  of  a  =  ±5.675°  was  found 
to  be  the  stall  angle  by  setting  the  2  fits  equal  to  each  other  and  solving  for  a. 
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Fig.  21  Canard  component  fits  in  local  blade  coordinate  system  for  Configuration  2,  Mach  0.65. 


Table  9  Static  aerodynamic  coefficient  fits  for  Canard  1  for  Configuration  2,  Mach  0.65 


a  range 

Cx0 

CXal 

CX„2 

Cxa3 

C*«4 

C  m 

7V  a 

CN„3 

r 

^rna  i 

r 

^rna  3 

Cla  3 

a  <  -5.675° 

-0.00911 

0.0 

1.399 

0.0 

-0.00130 

0.0 

-0.0160 

0.0 

-0.0396 

0.0 

-5.675°  <  a  <  5.675° 

0.00430 

0.0 

-1.416 

0.0 

64.751 

1.349 

-23.367 

1.205 

-20.410 

-1.055 

16.593 

a  >  5.675° 

-  -0.0443 

0.0 

-0.344 

0.0 

0.180 

0.0 

0.139 

0.0 

-0.223 

0.0 

Table  10  Static  aerodynamic  coefficient  fits  for  Canard  2  for  Configuration,  Mach  0.65 


a  range 

CXo 

Cxal 

CX„2 

CX„3 

C*«4 

C  m 

w  a 

CNU3 

r 

^rna  i 

r 

3 

Cla3 

a  <  -5.675° 

0.0440 

0.0 

0.378 

0.0 

0.191 

0.0 

0.149 

0.0 

-0.232 

0.0 

-5.675°  <  a  <  5.675° 

0.00430 

0.0 

-1.417 

0.0 

64.802 

1.349 

- 

1.204 

-20.315 

-1.055 

16.522 

a  >  5.675° 

0.00396 

0.0 

-1.297 

0.0 

-0.0110 

- 

-0.238 

0.0 

-0.0321 

0.0 

Table  1 1  Static  aerodynamic  coefficient  fits  for  Canard  3  for  Configuration  2,  Mach  0.65 


a  range 

cx0 

C*al 

CX„2 

C^a3 

C*a4 

C  m 

7V  a 

CNa3 

r 

^rnal 

r 

^rna3 

Cla  3 

a  <  -5.675° 

-0.00886 

0.0 

1.395 

0.0 

-0.00126 

0.0 

-0.0159 

0.0 

-0.0396 

0.0 

-5.675°  <  a  <  5.675° 

0.00430 

0.0 

-1.416 

0.0 

64.760 

1.349 

-23.374 

1.205 

-20.416 

-1.055 

16.597 

a  >  5.675° 

-0.0441 

0.0 

-0.349 

0.0 

0.180 

0.0 

0.139 

0.0 

-0.223 

0.0 

Table  12  Static  aerodynamic  coefficient  fits  for  Canard  4  for  Configuration  2,  Mach  0.65 


a  range 

CX0 

C^al 

CX„2 

Cxa3 

C^a4 

C  m 

7V  a 

CNa  3 

r 

^rnal 

r 

Cla  4 

cia3 

a  <  -5.675° 

0.0440 

0.0 

0.378 

0.0 

0.191 

0.0 

0.148 

0.0 

-0.232 

0.0 

-5.675°  <  a  <  5.675° 

0.00430 

0.0 

-1.417 

0.0 

64.791 

1.349 

-23.252 

1.204 

-20.311 

-1.055 

16.520 

a  >  5.675° 

0.00385 

0.0 

-1.295 

0.0 

-0.011 

0.0 

-0.0239 

0.0 

-0.0321 

0.0 

Having  completed  the  component  fits,  the  stability  of  the  projectile  and  the  canard  stall  angle  can 
be  investigated.  With  no  canard  deflection  (Fig.  22),  the  total  force  and  pitching  moment  should, 
and  do,  match  the  total  values  as  determined  by  CFD.  The  center-of  pressure  location  is 
predicted  to  be  just  under  0.5  calibers  behind  the  Xcg  for  angles  of  attack,  indicating  the  current 
design  is  indeed  stable.  Additionally,  canard  stall  occurs  at  approximately  a  =  ±8°  (projectile 
body  angle),  where  the  magnitude  of  the  canard  contribution  to  these  coefficients  becomes 
approximately  constant. 
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Fig.  22  Normal  force,  pitching  moment,  and  center-of-pressure  location  build  ups  using  CFD  component  fits  with 
0°  canard  deflection  for  Configuration  2,  Mach  0.65 

In  addition,  the  accuracy  of  assuming  that  the  interaction  effects  of  the  canards  on  the  fins  scales 
with  the  fins  (i.e.,  lumping  the  interaction  effects  into  the  fin  component),  as  was  done  in  the  first 
phase,  can  be  determined.  The  aerodynamic  coefficients  determined  from  the  CFD  simulation  of 
Configuration  2  (Fig.  22)  are  compared  to  those  predicted  when  the  scaling  factors  representative 
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of  Configuration  2  (fF  =  2.0 ,fc  =  0.7,xc  =  0.88,  Xcg  =  2.63)  were  applied  to  Configuration  la 
aerodynamic  data  (Fig.  23).  The  normal  force  and  pitching  moment  of  the  body  component 
(Fig.  24)  vary  between  the  CFD  solution  and  the  predictions  by  less  than  10%.  The  fin 
component  (Fig.  25)  predictions  under  predict  the  CFD  solution  by  25%  for  CN  and  15%  for  Cm. 
These  significantly  smaller  variations  indicate  that  combining  the  interaction  effects  with  the 
contribution  is  a  viable  option  during  the  initial  design  process. 


Projectile  Body  Angle-of-Attack  (deg) 


Fig.  23  Component  build  up  of  normal  force  and  pitching  moment  using  scaled  component  contributions 

ifF  =  2.0,  fc  =  0.7,  xc  =  0.88,  Xcg  =  2.63 )  from  CFD  analysis  of  Configuration  la  to  obtain  theoretical  values 
for  Configuration  2,  Mach  0.65 


Scaling  the  interaction  effects  caused  the  fin  effectiveness  to  be  under  predicted  whereas 
neglecting  it  caused  it  to  be  over  predicted.  The  under  prediction  is  preferable  as  it  is  easy  to 
make  a  projectile  less  stable.  Being  able  to  scale  the  interaction  effects  will  limit  the  number  of 
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simulations  and  grids  that  are  required  during  the  design  process  as  determination  of  interaction 
effects  would  require  additional  simulations  on  body-alone,  body-fin,  and  body-canard 
geometries.  The  differences  in  the  normal  force  and  pitching  moment  of  the  canard  component 
(Fig.  26)  appear  to  be  due  entirely  to  the  occurrence  of  canard  stall  in  Configuration  2.  Perhaps 
the  presence  of  the  additional  fins  (8  on  Configuration  2  vs.  4  on  Configurations  la  and  lb)  at 
subsonic  velocities  caused  feedback  upstream  to  the  canards. 


Fig.  24  Comparison  of  body  component  for  Configuration  2  as  determined  from  component  buildup  using 
Configuration  la  aerodynamics  with  form  factors  and  CFD  simulations,  Mach  0.65 


Fig.  25  Comparison  of  fin  component  for  Configuration  2  as  determined  from  component  buildup  using 
Configuration  la  aerodynamics  with  form  factors  and  CFD  simulations,  Mach  0.65 
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Fig.  26  Comparison  of  canard  component  for  Configuration  2  as  determined  from  component  buildup  using 
Configuration  la  aerodynamics  with  form  factors  and  CFD  simulations,  Mach  0.65 

Knowing  the  effect  of  canard  deflection  on  the  stability  of  the  projectile  was  also  desired.  Instead 
of  completing  additional  simulations  with  canard  deflection,  the  effect  of  canard  deflection  was 
once  again  estimated  by  using  total  canard  angle  of  attack  (8+  a)  to  determine  the  canard 
component  contribution  to  the  aerodynamics  of  the  airframe.  When  a  canard  deflection  angle  of 
8  =  4°  is  investigated  (Fig.  27),  canard  stall  is  predicted  to  occur  at  approximately  a  =  3°-4° 
body  angle  as  would  be  expected.  As  such,  total  Cm  doesn’t  continue  to  decrease  as  rapidly,  nor 
does  total  CN  increase  as  rapidly.  Examining  the  Cm  and  center-of-pressure  plots  in  Fig.  27,  the 
Configuration  2  airframe  reaches  a  trimmed  state  at  approximately  a  =  3°.  The  presence  of 
canard  stall  appears  to  cause  a  maximum  trim  angle  of  a  =  3°  even  as  the  deflection  angle 
further  increases  to  8=  6°  (Fig.  28).  The  actual  variation  of  canard  stall  angle  with  canard 
deflection  will  need  to  be  investigated  in  the  future  to  determine  what  the  maximum  trim  angle 
could  be. 
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Fig.  27  Normal  force,  pitching  moment,  and  center-of-pressure  location  build  ups  using  CFD  component  fits  with 
4°  canard  deflection  for  Configuration  2,  Mach  0.65 
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Fig.  28  Normal  force,  pitching  moment,  and  center-of-pressure  location  build  ups  using  CFD  component  fits  with 
6°  canard  deflection  for  Configuration  2,  Mach  0.65 

As  Configuration  2  was  determined  to  be  stable  for  small  to  moderate  canard  deflection  angles,  it 
is  likely  to  be  the  platform  of  choice  as  the  program  moves  forward.  Therefore,  additional 
knowledge  of  the  design,  such  as  equilibrium  roll  rate,  was  desired  so  that  necessary 
modifications  to  the  design  could  be  made  prior  to  moving  forward.  The  equilibrium  roll  rate  of 
the  current  design  is  a  balance  of  the  roll  torque  being  generated  by  the  bend  of  the  fin  tips  and 
the  roll  damping  generated  by  the  body  as  well  as  the  size  of  the  fins  and  canards.  The  roll  torque 
coefficient  at  a  =  0°  was  determined  from  the  static  computation.  However,  a  separate 
simulation,  a  rotating  reference  frame  computation  using  the  technique  in  DeSpirito  et  al,15  was 
required  to  determine  the  roll  damping  coefficient.  From  this  computation,  the  roll  damping 
coefficient,  Cip,  for  the  complete  configuration  (body-tail-canards)  was  determined  to  be  -10.392. 
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Equation  28  was  used  to  predict  the  equilibrium  roll  rate,  peq,  of  6.34  Hz.  While  this  roll  rate  is 
about  double  the  desired  roll  rate,  it  remains  a  good  place  to  start. 


Peq  — 


Cj_2V 
C,  D 

Lp 


(28) 


While  the  primary  purpose  of  this  investigation  was  to  design  a  stable  platform  for  the  LCHAW 
moving  target  mission  program,  another  research  goal  was  to  begin  to  better  understand  the 
interaction  of  the  canard  tip  vortices  with  the  fins.  Figure  29  shows  how  the  streamlines  that 
begin  at  the  canards  with  5  =  0°  change  as  the  angle  of  attack  increases  from  a  =  0°  to  a=  14°. 
For  5  =  0°  and  a  =  0°,  there  is  a  complete  lack  of  canard  tip  vortices  and  interference  with  the 
fins.  The  lack  of  tip  vortices  and  interference  continues  at  small  angles  of  attack.  As  a  increases 
to  3°,  canard  tip  vortices  are  starting  to  develop.  It  is  not  until  a  increases  beyond  5°  that  the 
vortices  appear  to  twist  enough  to  interfere  with  the  fins.  Interference  is  minimal  at  a=  5°,  but  is 
occurring  by  a  =  8°.  Figure  30  shows  the  a  =  10°  case  from  the  rear  of  the  projectile. 

The  necessity  of  accounting  for  the  interaction  effects  for  all  4  canards  becomes  quickly 
apparent.  Obviously,  not  all  of  the  fins  are  affected,  especially  those  fins  on  the  windward  side  of 
the  projectile.  However,  only  some  of  the  fins  are  affected  by  the  tip  vortices;  the  tip  vortices  of 
the  top  canards  are  swept  away  from  the  fins  when  the  projectile  is  pitched  up.  The  streamlines 
from  inboard  of  all  4  canards  appear  to  affect  both  the  body  and  the  fins.  The  exact  effects  will 
have  to  be  investigated  further,  especially  the  effects  when  the  canards  are  deflected. 
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Fig.  29  Pitch  plane  Mach  number  contours  with  pressure  coefficient  contours  on  fins  and  canards  (S  =  0 °)  with 
streamlines  at  various  angles  of  attack  for  Configuration  2,  Mach  0.65 


42 


Fig.  30  Individual  canard  streamlines,  rear  view,  for  Configuration  2,  Mach  0.65 

4.2.2  Parameter  Sensitivity  Study 

Because  there  were  differences  between  the  predicted  and  actual  aerodynamic  performance  of 
the  round,  it  was  decided  to  use  the  parameter  sensitivity  routine  to  determine  if  any  further 
modifications  to  the  geometry  could  provide  additional  maneuverability  through  an  increase  to 
the  lift-to  drag  ratio.  Thus,  the  aerodynamic  data  for  Configuration  2  were  used  as  the  input  for 
the  optimization  and  the  parameters  of  interest  (fin  and  canard  area  through  fF  and  fc )  varied. 
Unfortunately,  canard  axis  location  was  not  shifted  as  the  canard  area  (i.e.,  span  as  the  chord  is 
fixed  due  to  geometric  constraints)  was  varied,  which  may  not  provide  for  ideal  results.  These 
results  are  provided  in  Fig.  31  for  Configuration  2  with  a  reduction  in  fin  area  of  20%  (fF  =  0.8), 
an  increase  in  canard  area  of  50%  ( fc  =  1.5),  and  5  =  4°.  This  specific  configuration,  which  is 
likely  marginally  stable,  provides  a  maximum  lift-to-drag  ratio  around  3,  and  produces  a  trim 
angle  of  approximately  6°.  The  trends  with  fin  size,  canard  size,  and  canard  deflection  are 
similar  to  those  for  Configurations  la  and  lb.  These  data  indicate  that  some  additional  fin  area 
reduction  and  canard  area  addition  may  enhance  maneuverability  (i.e.,  more  normal  force,  larger 
trim  angle,  etc.).  However,  some  margin  is  warranted  to  account  for  uncertainty.  Thus, 
Configuration  2  is  chosen  as  the  airframe  for  further  consideration. 
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Fig.  31  Optimization  analysis  using  Configuration  2  aerodynamic  data  with  /V  =  0. 8,  f  =  1.5,  and  5  =  4° 

4.2.3  Linear  Flight  Mechanics  and  Control 

A  detailed  study  of  the  controlled  flight  performance  was  undertaken  based  on  the  favorable 
aerodynamic  characteristics  of  Configuration  2.  The  linear  flight  mechanics  were  implemented  in 
simulation  and  driven  with  the  appropriate  aerodynamic,  mass  property,  and  actuator  data.  The 
flight  control  commands  canard  deflection  to  achieve  a  desired  lateral  acceleration  while 
minimizing  angular  rate. 

Some  linear  simulation  results  in  the  pitch  plane  for  Configuration  2  are  provided  in  Fig.  32.  The 
top  left  plot  shows  that  the  control  sufficiently  reduces  the  pitch  rate  to  the  reference  pitch  rate 
(i qr )  of  zero.  The  measured  pitch  rate  (qM)  and  actual  pitch  rate  (q)  are  equivalent  since  perfect 
feedback  was  used  in  the  simulations.  A  reference  pitch  acceleration  (wr)  of  40m/s2  was  desired 
and  the  response  (measured  pitch  acceleration  (wM)  and  actual  pitch  acceleration  (w)  equivalent 
since  perfect  feedback)  suggests  good  performance  since  the  rise  time  is  less  than  0.1s.  The 
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commanded  and  actual  pitch  deflections  ( 8qc  and  Sq,  respectively)  to  achieve  these  pitch  rate 
and  pitch  acceleration  responses  are  also  shown  in  Fig.  32.  Approximately  9° of  deflection 
yielded  40m/s2  of  acceleration. 


Fig.  32  Response  and  commands  in  pitch  from  linear  modeling  and  simulation 

4.2.4  Nonlinear  Flight  Mechanics  and  Control 

Simulations  of  the  controlled  flight  were  performed  on  Configuration  2  using  the  nonlinear  flight 
mechanics  models.  A  target  was  placed  approximately  1000m  downrange  moving  about  5m/s  in 
the  crossrange  direction.  Uncertainties  or  variations  were  placed  on  the  target  characteristics, 
launch,  and  projectile  technologies  to  account  for  real-world  effects.  Uncompensated  launch 
error  sources  (e.g.,  muzzle  velocity  variation,  gun  point  angle,  and  initial  angular  disturbances) 
were  modeled.  The  roll  angle,  angular  rate  sensor,  accelerometer,  and  line-of-sight  rate  feedback 
were  corrupted.  Nominal  actuator  characteristics,  projectile  mass  properties,  and  aerodynamics 
were  perturbed.  Monte  Carlo  techniques  were  applied  in  the  simulations. 
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The  trajectory  and  aerodynamic  angle  histories  for  a  representative  flight  are  presented  in  Fig.  33. 
The  guided  projectile  is  the  solid  blue  trace  and  the  target  is  the  dashed  red  line  in  the  trajectory 
plots.  The  projectile  moves  about  30  m  off  the  line-of-fire  in  the  crossrange  and  altitude 
directions  due  to  maneuvers  necessary  to  intercept  the  target  and  the  action  of  gravity.  These  data 
indicate  a  successful  moving  target  intercept  with  small  point-of-closest-approach.  The 
maneuvers  required  for  this  engagement  are  well  within  the  bounds  of  the  Configuration  2 
airframe  since  the  pitch  and  yaw  angles-of-attack  are  less  than  1°  throughout  this  flight. 


Fig.  33  Trajectory  and  aerodynamic  angle  histories  for  representative  flight  of  Configuration  2 

The  objective  of  the  flight  control  is  to  track  specific  roll  angles,  null  out  angular  rates,  and 
achieve  desired  lateral  accelerations.  The  performance  of  the  controller  for  Configuration  2  using 
the  nonlinear  models  for  this  representative  flight  is  provided  in  Fig.  34.  Each  plot  features  the 
true  value,  the  measured  (corrupted  truth)  value,  and,  for  all  plots  except  roll  angle,  the  desired 
value.  The  roll  angle  is  controlled  to  within  10°  of  the  reference  point.  Inspection  of  all  angular 
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rate  histories  (roll,  pitch,  and  yaw)  clearly  shows  good  performance  because  rates  are  near  zero 
for  the  majority  of  the  flight.  Pitch  and  yaw  acceleration  tracking  is  satisfactory.  Measurement 
error  is  evident  in  all  plots,  although  specifically  called  out  for  the  roll  angle. 
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Fig.  34  Controlled  response  histories  in  roll,  pitch,  and  yaw  for  representative  flight  of  Configuration  2 
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The  individual  canard  histories  from  this  representative  flight  are  given  in  Fig.  35.  The 
commands  are  shown  along  with  the  true  value.  A  first-order  model  with  a  bias  causes  the 
command  and  actual  values  to  differ.  Some  action  is  evident  early  in  the  flight  to  initially  control 
the  roll  angle  and  angular  rates.  Only  small  variation  is  required  throughout  the  majority  of  the 
flight  to  meet  the  lateral  acceleration  requirements  for  intercepting  the  target. 
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Fig.  35  Canard  deflection  histories  for  representative  flight  of  Configuration  2 

Canard  deflection  saturation  (i.e.,  reaching  a  maximum  value,  for  example,  of  10°  in  Fig.  35)  in 
turn  saturates  the  lateral  acceleration  (i.e.,  constant  v  or  w  in  Fig.  34,  for  example,  of  50m/s2)  at 
the  limit  based  on  the  airframe  characteristics.  The  lack  of  saturation  in  the  results  for  the  canard 
commands  and  controlled  responses  in  Figs.  34  and  35  confirms  that  the  design  methodology 
yielded  an  airframe  capable  of  successful  intercept  for  this  class  of  launch  and  target  conditions. 


48 


5.  Conclusions 


The  primary  purpose  of  this  investigation  was  to  design  a  stable  platform  for  the  LCHAW 
moving  target  mission  program.  This  study  demonstrated  the  need  to  include  canard  interaction 
effects  when  determining  the  stability  a  of  body-fin-canard  projectile  configuration.  The  body- 
fin-canard  projectile  that  had  been  designed  using  the  aerodynamic  data  obtained  for  the  body- 
fin  portion  of  the  projectile  without  canards  and  then  adding  the  canard  data  from  a  separate 
source  was  found  to  be  statically  unstable  across  the  Mach  number  regime  of  interest  when  the 
aerodynamic  data  for  the  components  were  obtained  from  a  single  CFD  simulation.  The 
component  aerodynamic  data  from  the  single  CFD  simulation  (i.e.,  canard  interaction  effects 
included)  was  used  in  a  parameter  sensitivity  study  to  determine  a  new  configuration  that  would 
meet  stability  requirements.  The  results  of  the  parameter  sensitivity  study  plus  geometric 
constraints  were  used  to  design  a  second  configuration  for  consideration  as  the  demonstration 
platform  for  the  LCHAW  moving  target  mission  program. 

Using  form  factors  on  data  that  already  included  the  canard  interaction  effects  produced 
reasonable  results — the  updated  configuration  remained  statically  stable  with  a  trim  angle  of 
approximately  4°  with  deflected  canards.  Comparing  the  computational  aerodynamic  coefficients 
to  the  predicted  values  proved  that  incorporating  interaction  effects  into  the  component 
aerodynamic  coefficient,  even  if  scaling  was  necessary,  was  a  reasonable  assumption  during  the 
design  process.  A  complete  aerodynamic  model  was  developed  for  this  second  configuration, 
including  roll  damping  coefficient,  at  a  single  Mach  number  using  CFD  simulations.  The 
aerodynamic  model  was  used  as  input  for  both  linear  and  non-linear  flight  mechanics  and  control 
simulations.  These  simulations  showed  that  this  configuration  would  meet  maneuverability  and 
range  requirements  as  well  as  stability  requirements. 

During  the  course  of  the  design  process,  a  number  of  interesting  phenomena  were  discovered 
that  need  to  be  investigated  further.  A  most  basic  concern  was  the  inability  of  the  ballistic  fits  to 
predict  C*  and  Clo  accurately.  Additional  work  on  the  fits  should  clear  up  this  issue.  Currently, 

only  configurations  with  non-deflected  canards  were  considered.  However,  there  are  2  items  that 
lend  themselves  to  a  more  in  depth  investigation  once  canard  deflections  are  considered.  First, 
the  effect  of  canard  deflection  on  canard  stall  angle  and  projectile  trim  angles  will  need  to  be 
investigated.  Second,  the  interaction  effects  between  the  canard  and  both  the  body  and  the  fins 
when  the  canards  are  deflected  will  need  further  research.  It  is  likely  that  tracing  the  canard 
streamlines  emanating  from  both  the  inboard  and  tip  locations  will  be  necessary  to  determine  the 
full  extent  of  the  interactions.  The  interaction  effects  will  likely  be  critical  in  understanding  and 
predicting  the  projectile  maneuver  behavior,  especially  at  higher  angles  of  attack. 
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Appendix  A.  Transformations 
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Transformations  between  coordinate  systems  were  an  integral  part  of  the  analysis  presented  in 
the  main  body  of  this  report.  The  transformations  ensured  that  the  aerodynamic  coefficient 
database  created  from  the  computational  fluid  dynamics  (CFD)  analysis  was  properly  utilized  in 
both  the  optimization  routine  as  well  as  in  the  flight  mechanics  and  control  analysis. 

The  first  transformation  was  completed  to  move  from  the  CFD  coordinate  system  to  the 
projectile  body  coordinate  system.  The  CFD  coordinate  system  is  the  right-hand  coordinate 
system  with  the  x-axis  aligned  with  the  centerline  of  the  projectile  positive  from  nose  to  tail  and 
the  z-axis  positive  up,  as  defined  in  Fig.  A-l.  The  body-fixed  coordinate  system  is  defined  in 
Fig.  A-2.  A  major  difference  in  these  2  coordinate  systems  stems  from  the  definition  of  the 
velocity  vector;  in  the  CFD  coordinate  system  the  air  moves  over  the  projectile,  while  the 
projectile  moves  through  the  air  in  the  body-fixed  coordinate  system.  Thus,  VCG /t  =  —Vwind. 


-CFD 


Fig.  A-l  CFD  coordinate  system  and  aerodynamic  angles 
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B 


CG/I 


Fig.  A-2  Body-fixed  coordinate  system  and  aerodynamic  angles 

The  transformation  matrix  to  convert  from  CFD  to  projectile  body  coordinates,  TCB,  is  given  by 
Eq.  A-l.  Equation  A-2  is  then  used  to  convert  between  the  2  coordinate  systems 


-l 

0 

CB 

~ 

0 
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.  0 

0 

X 

X 

Y 

=  tcb 

Y 

2. 

B 

2. 

0 

0 


(A-l) 


(A-2) 


where 


CFD 


are  the  3  components  of  the  force  or  moment  coefficients  or  velocity.  Equation  A-2 


only  transforms  the  velocity  vector  from  CFD  coordinates  to  body  coordinates;  the  velocity  vector 
remains  the  wind  velocity  vector.  To  complete  the  transformation  to  body  coordinates  as  defined 
in  Fig.  A-2,  Vwind  must  be  converted  to  VCG /,  using  Eq.  A-3. 


Vcg/i  — 


-10  0 
0-10 
0  0-1 


FWI- 


wind 


(A-3) 


It  was  determined  that  the  canard  components  should  be  analyzed  in  the  local  blade  coordinate 
system  rather  than  in  the  projectile  body  coordinate  system.  As  such,  the  force  and  moments 
coefficients  for  each  of  the  canards  must  further  be  transformed  from  the  projectile  body 
coordinate  to  the  local  blade  coordinate  system.  The  local  blade  coordinate  system  is  defined  as 
the  x-axis  parallel  to  xB.  The  y-axis  runs  along  the  span  of  the  canard  from  the  root  to  the  tip. 
The  z-axis  runs  perpendicular  to  follow  the  right  hand  rule.  The  transformation  matrix  to  convert 
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from  the  projectile  body  coordinate  system  to  the  local  blade  coordinate  system  of  the  ith  canard, 
Tbm.,  is  given  by  Eq.  A-4,  where  0M.  is  the  roll  orientation  of  the  given  canard.  The  roll 

orientation  of  the  canards  is  determined  by  looking  from  the  rear  of  the  projectile  (Fig.  A-3).  The 
roll  orientation  of  Canard  4  (i.e.  M4 )  is  <pM.  =  45°  for  this  study.  In  Fig.  A-3,  the  sign  convention 

for  canard  deflection  is  for  the  trailing  edge.  The  aerodynamic  coefficients  in  the  local  blade 
coordinate  system  of  the  ith  canard  are  then  determined  by  Eq.  A-5. 


Fig.  A-3  Definition  of  canard  orientation  when  looking  from  the 
rear  of  the  projectile  relative  to  body  coordinates 

The  final  transformations  completed  in  this  paper  are  to  change  the  aerodynamic  coefficients  in 
the  projectile  body  frame  to  ballistic  coefficients  (Fig.  A-4).  To  obtain  the  ballistic  coefficients, 
only  the  sign  on  vertical  force  (i.e.,  normal  force)  and  axial  force  must  be  inverted  as  CN  is 
defined  positive  up  and  CA  positive  aft  for  positive  a,  such  that 


CN  — 

~CzB 

(A-6) 

CA  = 

~CXB 

(A-7) 

56 


Within  this  paper,  CA  and  Cx  are  used  interchangeably  to  indicate  axial  force  coefficient  as 
Eq.  A-7  was  applied  early  in  the  transformation  process  without  consequence. 


Fig.  A-4  Definition  of  ballistic  coefficients  and  aerodynamic  angles  in  the  body  coordinate  system 
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Appendix  B.  Coefficient  Fits  for  Configurations  la  and  lb 
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Best  fits  were  determined  for  each  of  the  components  for  both  Configurations  la  and  lb  at  all  3 
Mach  numbers.  The  resulting  fits  were  checked  for  accuracy  against  the  CFD  for  both  the 
individual  component  and  for  the  sum  of  the  components.  The  value  of  the  total  fit  was  obtained 
by  summing  the  component  fits;  the  total  CFD  value  was  obtained  directly  from  the  CFD  (i.e., 
the  CFD  components  were  not  summed).  Using  the  total  CFD  value  allowed  for  an  additional 
check  that  all  component  values  were  being  correctly  utilized. 

Figures  B-l  and  B-2  show  the  resulting  fits  for  Configuration  la  at  Mach  0.5  for  the  body  and  fin 
components  and  each  of  the  canard  components,  respectively.  Figure  B-3  shows  the  results  of  the 
component  buildup. 


Fig.  B-l  Body  and  fin  component  fits  for  Configuration  la,  Mach  0.5 
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Fig.  B-2  Individual  canard  fits  in  local  canard  coordinate  system  for  Configuration  la,  Mach  0.5 
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Fig.  B-3  Component  fit  build  up  comparison  to  total  CFD  coefficients  for  Configuration  la,  Mach  0.5 


Figures  B-4  and  B-5  show  the  resulting  fits  for  Configuration  la  at  Mach  0.65  for  the  body  and 
fin  components  and  each  of  the  canard  components,  respectively.  Figure  B-6  shows  the  results  of 
the  component  buildup. 
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Fig.  B-4  Body  and  fin  component  fits  for  Configuration  la,  Mach  0.65 
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Fig.  B-5  Individual  canard  fits  in  local  canard  coordinate  system  for  Configuration  la,  Mach  0.65 
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Fig.  B-6  Component  fit  build  up  comparison  to  total  CFD  coefficients  for  Configuration  la,  Mach  0.65 


Figures  B-7  and  B-8  show  the  resulting  fits  for  Configuration  la  at  Mach  0.8  for  the  body  and  fin 
components  and  each  of  the  canard  components,  respectively.  Figure  B-9  shows  the  results  of  the 
component  buildup. 
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Fig.  B-7  Body  and  fin  component  fits  for  Configuration  la,  Mach  0.8 
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Fig.  B-8  Individual  canard  fits  in  local  canard  coordinate  system  for  Configuration  la,  Mach  0.8 
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Fig.  B-9  Component  fit  build  up  comparison  to  total  CFD  coefficients  for  Configuration  la,  Mach  0.8 


Figures  B-10  and  B-ll  show  the  resulting  fits  for  Configuration  lb  at  Mach  0.5  for  the  body  and 
fin  components  and  each  of  the  canard  components,  respectively.  Figure  B-12  shows  the  results 
of  the  component  buildup. 
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Fig.  B-10  Body  and  fin  component  fits  for  Configuration  lb,  Mach  0.5 
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Fig.  B-l  1  Individual  canard  fits  in  local  canard  coordinate  system  for  Configuration  lb,  Mach  0.5 
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Fig.  B-12  Component  fit  build  up  comparison  to  total  CFD  coefficients  for  Configuration  lb,  Mach  0.5 


Figures  B-13  and  B-14  show  the  resulting  fits  for  Configuration  lb  at  Mach  0.65  for  the  body 
and  fin  components  and  each  of  the  canard  components,  respectively.  Figure  B-15  shows  the 
results  of  the  component  buildup. 
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Fig.  B-13  Body  and  fin  component  fits  for  Configuration  lb,  Mach  0.65 
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Fig.  B-14 


Individual  canard  fits  in  local  canard  coordinate  system  for  Configuration  lb,  Mach  0.65 
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Fig.  B-15  Component  fit  build  up  comparison  to  total  CFD  coefficients  for  Configuration  lb,  Mach  0.65 


Figures  B-16  and  B-17  show  the  resulting  fits  for  Configuration  lb  at  Mach  0.8  for  the  body  and 
fin  components  and  each  of  the  canard  components,  respectively.  Figure  B-18  shows  the  results 
of  the  component  buildup. 
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Fig.  B-16  Body  and  fin  component  fits  for  Configuration  lb,  Mach  0.8 
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Individual  canard  fits  in  local  canard  coordinate  system  for  Configuration  lb,  Mach  0.8 
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Fig.  B-18  Component  fit  build  up  comparison  to  total  CFD  coefficients  for  Configuration  lb,  Mach  0.8 
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List  of  Symbols,  Abbreviations,  and  Acronyms 

wM  measured  pitch  acceleration  (m/s2) 

wr  reference  pitch  acceleration  (m/s2) 

Clp  roll  damping  moment  coefficient 

CA  axial  force  coefficient 

CN  normal  force  coefficient 

Ci  axial  moment  (roll  torque)  coefficient 

Cm  pitching  moment  coefficient 

qM  measured  pitch  rate  (m/s) 

qr  reference  pitch  rate  (m/s) 

w  actual  pitch  acceleration  (m/s2) 

xCref  canard  hinge  location  at  which  canard  aerodynamic  coefficients  were  determined 

(cal.) 

B  body 

C  canard 

Cp  pressure  coefficient 

D  reference  diameter  (m) 

dy  first  cell  spacing  off  wall  (m) 

F  fin 

fc  canard-scaling  factor 

fp  fin-scaling  factor 

2 

Ix  axial  moment  of  inertia  (kg-m  ) 

Iy  transverse  moment  of  inertia  (kg-m  ) 

k  turbulent  kinetic  energy 

L/D  lift-to-drag 

M  Mach  number 
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m 

Peq 

q 

R 

s 

V 

xc 

xcg 

y+ 

a 

(%trim 

5 

£ 

P 

1- D 

2- D 

3- D 

ARDEC 

ARL 

CAS 

CFD 

CFL 

DOD 

DSRC 

HPCMP 

LCHAW 


mass  (kg) 

equilibrium  spin  rate  (Hz) 

i  „ 

dynamic  pressure,  -  pV  ,  actual  pitch  rate  in  flight  dynamic  simulations 
undamped  eddy  viscosity 
aerodynamic  reference  area, 

velocity  magnitude  (m) 

canard  hinge  location  forward  of  Xcg  (cal.) 

center  of  gravity  location  (m) 

non-dimensional  wall  spacing 

angle  of  attack  (degrees) 

trim  angle  of  attack  (degrees) 

canard  deflection  angle  (or  sin  a  where  defined  in  total  coefficient  equations) 
dissipation  rate 
density  (kg/m3) 

1 - dimensional 

2- dimensional 

3- dimensional 

Armament  Research,  Development,  and  Engineering  Center 

Army  Research  Eaboratory 

control  actuation  system 

computational  fluid  dynamics 

Courant-Freidrich-Eewy 

Department  of  Defense 

DOD  Supercomputing  Resource  Center 

High-Performance  Computing  Modernization  Program 

Low  Cost  Hyper  Accurate  Weapon 
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MPI 

PVI 

RANS 

YAPP 


Message  Passing  Interface 
Parallel  Virtual  Machine 
Reynolds-averaged  Navier-Stokes 
Very  Affordable  Precision  Projectile 
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